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Abstract 

We find an upper bound to the maximal number of hmit cycles, 
which bifurcate from a hamiltonian two-saddle loop of an analytic 
vector field, under an analytic deformation. 
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1 Introduction 

Consider a A^-parameter analytic family of analytic plane vector fields Xx, 
X G (M^, 0), such that Xq has a A;-saddle cycle (a hyperbolic fc-graphic) F^, as 
on figjT] The cychcity Cycl{Tk, Xx) of Fjt is, roughly speaking, the maximal 
number of limit cycles of Xx which tend to as A 0. The first results on 
the cychcity of one-saddle connection (also called homoclinic connection) go 
back to Andronov and Leontovich in 1937 (but they were pubhshed only in 
1959 P). The cychcity Cycl{Ti) has been studied later in full generality by 
Roussarie [El HI], see also [HI [12] for an extensive list of references. 

The main technical tool of the Roussarie's method is an asymptotic ex- 
pansion of the Dulac map (transport map near the saddle point) 

X dx{x) 

in terms of x^ ,x^uj{x,e) where uj{x,e) is the so called Ecalle-Roussarie com- 
pensator 

x^"^ — 1 

uj{x,e) = ,uj{x,0) = - lnx,x e (M+,0) 
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and e = e{\) is the trace of the vector field Xx at the saddle point. Let 
Pa be the Poincare first return map, associated to Fi and Xx. The most 
delicate case to be studied is when Fi is of infinite co-dimension (Pq = id). 
For A ~ the map Pa is composed by a Dulac map (near the saddle point) 
and an analytic map (the transport map along the homoclinic orbit). The 
usual derivation-division algorithm then provides an upper bound for the 
cyclicity in terms of the number of the coefficients of the asymptotic series of 
the displacement map Pa — id, which vanish as A = 0. The same method was 
applied more recently to one-parameter deformations of Hamiltonian two- 
saddle loops (called also heterochnic Hamiltonian connections), under the 
non-generic assumption that one of the separatrices of F2 remains unbroken 
[His]. Recall that a fc-saddle cycle F^ is said to be Hamiltonian, provided that 
there is a neighborhood of F^ in which Xq allows an analytic first integral 
with only Morse critical points. 

The purpose of the present paper is to extend these results to the case of 
an arbitrary analytic perturbation of a Hamiltonian two-loop V2, having two 
hyperbolic equilibrium points. 

Our approach is different, as we do not use the asymptotic series of the 
corresponding Dulac maps d\ ^^'^ d\, shown on figji]). Recall that in the 
one-parameter case A = e G (M, 0), the displacement function d]. — can be 
approximated by an appropriate Abelian integral /(.) (or more generally, an 
iterated path integral) depending on a parameter t as follows 

d^(t)-d2(t), = e'^/(t) + ...,£ ~0. 

Therefore to count the zeros of dl — d^ (corresponding to limit cycles) it is 
enough to count the zeros of /(.), which can be done by making use of the 
so called "Petrov trick" (based on the argument principle), see [U] and the 
references given there. 

The above considerations hold true at least far from the singular points 
of the vector field Xa. As discovered in [H [3], however, not all limit cycles 
in a neighborhood of a two-saddle loop can be approximated in such a way. 
The missing "alien" limit cycles are moreover non-avoidable in generic N- 
parameter deformations Xa with X > 4. For this reason, we apply the 
argument principle directly to the displacement map dl{.) — d'^{.)s in an 
appropriate complex domain, in order to obtain an estimate to the number 
of its complex zeros (corresponding to complex limit cycles). 

The main technical result of the paper is Lemma [2| which claims that the 
zero locus of the imaginary part of the Dulac map is a real analytic curve of 
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Figure 1: Hamiltonian /c-saddle cycles 



= C at the origin. This makes possible to investigate the number of the 
zeros of the imaginary part of d\ — d\ along the zero locus of the imaginary 
part of d\. Indeed, the intersection numbers of two analytic curves is easily 
computed. The proof of our main result, Theorem [TJ is then completed by 
making use of the Petrov trick. 

In the course of the proof of Theorem [T] we assume, for the sake of sim- 
plicity, that our deformations depend on a single small parameter e. General 
multi-parameter deformations A Xx of Xq are then studied along the same 
lines, as it follows from the Hironaka's desingularization theorem. We explain 
this in Appendix B.2, see Theorem [4j 

Deformations of an arbitrary (possibly non-Hamiltonian) two-loop of in- 
finite co-dimension (Pa = id) can be studied in a similar way, and will be 
considered in another paper. 



2 Description of the result 

Let Xq be a real plane vector field. Recall that a polycycle of Xq is a topolog- 
ical polygon composed of separatrices and singular points. A fc-saddle cycle 
of Xq (or a hyperbohc fc-graphic) denoted F^, is a polycycle composed of k 
distinct saddle-type singular points pi,p2, ■ ■ ■ ,Pk, Pk+i = Pi and separatrices 
(heteroclinic orbits) connecting pi to Pi+i. Let cr be a segment transversal 
to the polycycle. The /c-saddle cycle is said to be Hamiltonian, provided 
that Xq has an analytic first integral / having Morse critical points at pi. It 
follows that Ffc bounds an annulus of periodic orbits {{x,y) : f{x,y) = t}t 
of Xq. Thus, a Hamiltonian 0-saddle cycle is simply a center, a Hamiltonian 

1- saddle cycle is a homoclinic loop bounding a period annulus, a Hamiltonian 

2- saddle cycle is a double heterochnic loop bounding a period annulus etc., 
see figjl} 
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One can find a "tubular neighborhood" U C = of 



rfcC{(a;,|/)eM2:/(x,y) = 0} 

such that 

• U is compact smooth manifold with a (real three-dimensional) border. 

• / is analytic in some neighborhood of U 

• the border dlJ is transversal to the complex curves {{x,y) e : 
f{x,y) = t}, provided that |t| is sufficiently small. 

• the intersection of U with the singular fiber {{x,y) E : f(x,y) = 0} 
is a union of k Riemann surfaces Di, each of them homeomorphic to 
an open disc. intersects transversally -Dj+i at p,, i — 1 . . .k, and 
DiDDj = for ^ 1. 

It follows that 

f:U^C 

defines a locally trivial fibration over a punctured neighborhood of the origin 
in C, and each fiber 

F, = Un{{x,y)eC^: f{x, y)=t},tj^O (1) 

is homeomorphic to a genus one surface with k punctures. 

A one-parameter analytic deformation of Xq is a a family of real- 
analytic plane vector fields, depending analytically on a real parameter e e 
(R, 0), and defined in a suitable neighborhood of the /c-saddle cycle r^. The 
corresponding fohation J^^ has an extension in a complex domain denoted by 
the same letter, and defined by 

df + euj,^0 (2) 

where ujs = P{x,y,e)dx + Q{x,y,e)dy is a one-form, and P,Q are real- 
analytic in x,y,£ in a neighborhood of F^. 

Parameterize the segment a by the "synchronized" local variable t = /|cr 
and let ^(t) C > be the continuous family of periodic orbits of Xq 
which tend to the poly cycle Tk C Fq as t tends to 0. To the family {'y{t)} 
we associate the trivial first return map 

Pq : a ^ a, Pq — id 



5 



which allows an analytic continuation for e 7^ to a first return map 



t ^ P,{t) = t + e'^Md{t) + ... 

The dots above mean a function in t, e which, for every fixed t such that -Pe(t), 
£ ~ is defined, is of the type 0(e'^"'"^). The so called Poincare-Pontryagin 
function function may be exphcitly computed, see [TBI [31 [B]. 

More generally, let 7(t) C Ft be any continuous family of closed loops 
intersecting the cross-section a. For \e\ sufficiently small we define in a similar 
way the holonomy map 

hl^\ (J ^ a 

related to the family of loops {7(t)}t and the deformed fohation JF^. By 
analogy to the Poincare return map we have 

=t + £<^M,(t) + ... 

t ^ hf^it) =t + e''Mdit) + ... 

where the meaning of the dots is as before, and the number d depends on 
{7(t)}t and . 

The holonomy map /i^ depends on the choice of a. In contrast to this, 
the Poincare-Pontryagin function does not depend on the cross-section a, 
it depends on the free homotopy class of the loop 7(t) only. Further, it can 
be expressed in terms of iterated path integrals of length at most d, along 
suitable meromorphic differential one-forms. It satisfies therefore a linear 
differential equation which has a Fuchs type singularity at t = 0, see [71 [9]. 
Thus, the leading term of has the form 

where p is an eigenvalue of the indicial equation of the Fuchsian equation 
related to the regular singular point t = 0. 



Definition 1 We shall callp the characteristic number of the holonomy map 

^{hf^) =P- 



hl^ and denote 
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In the Hamiltonian case the number //(P^) is rational, because the corre- 
sponding monodromy operator is quasi-unipotent (TOj. 

To formulate the main result of the paper consider, more specifically, the 
case k = 2 { a double heterochnic loop). As it follows from [9], the function 
Md is in fact an Abelian integral and can be written as 



for suitable analytic one-form a). 

Let (t) be two continuous famihes of closed loops vanishing at the 

saddle points pi andp2- We suppose that orientations of the loops "agree" in 
the sense that the intersection indices of the homotopy classes of 6i{t),62{t) 
with the homotopy class of the periodic orbit ■yit) is one and the same. 

The cychcity Cycl{T2,^e) of the 2-saddle cycle T2 with respect to the 
deformed foliation Ts is the maximal number of limit cycles which bifurcate 
from r2 near 6 = 0, see pjj for a precise definition. An upper bound for the 
cyclicity Cyd(T2, J-'e) is given in terms of the characteristic numbers of the 
holonomies associated to T2 as follows 

Theorem 1 



Cyd{T2, -F,) < 1 + u{P,) + max{z/(/if J, z/(/if J} + u{hl o h^J. (5) 



It is tempting to conjecture that in general the cyclicity Cycl(T2,J^£) is 
bounded by a similar expression in terms of the characteristic numbers of 
the holonomies associated to Ffc. Indeed, in the homoclinic case, = 1, by 
repeating the proof of Theorem[T] one obtains 




Cyd{T^,J^,)<u{P,) + u{hs,). 



We have, typically 



Psit) = t + eMi{t) + ... 



and if 



Mi(t) = hit) log(t) + f2{t) ^ 0, = 0{tn, f2{t) = 0{n 



for some analytic functions /i,/2, then 



Cyd(Ti, Te) < min{p, q) + q. 
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By the Roussarie's theorem [T71 Theorem C] the exact upper bound in a 
real domain in this case is 2p if p < g, and 2g — 1 if p > g. This suggests 
that the bound of Theorem [T] can be improved. In fact, the bound ^ holds 
true for the number of complex limit cycles accumulating on r2 in a suitable 
neighborhood of it. 

Example. Suppose that d = 1 in (3), that is to say Mi(t) = Xy(() ^^o where 
{7(t)}t>o is the family of real perioaic orbits of JFq = {df = 0}. Then we 
have 

Mi(t) = (A(t) + /2(t))log(t) + /3(t) 

where the functions /i, /2, /s are analytic in a neighborhood of t = 0, 



fl{t) = / UJ0,f2{t) = UJO 
JSlit) JS2{t) 

hs,{t) = t + ehit) + ...,hs,{t)=t + £f2{t) + ... 
hs, o hs,{t) = t + e{Mt) + Mt)) + . . . 
Suppose further that 

Mt) = 0{tn, hit) = oit^^), hit) = o{r^), m + hit) = oit"). 

Theorem [T] implies that the cyclicity of r2 is bounded by 

1 + min{p, q} + max{pi,p2} + q- 

In the case pi = p2 = q = p iov instance, this gives 

Cyd(T2,J^e) <l + 3p. 

In this situation, and under the strong hypothesis that one of the connections 
of r2 remains unbroken, it has been proved in [U, Theorem 8] that 

Cyd{T„J^,)<2p-l + ^^^^. 



3 The Dulac map 

Let be a real analytic foliation defined by ([2]) in a neighborhood of a 
hyperbohc Morse critical point of the function /. For all sufficiently small 
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Figure 2: The Dulac map 



\e\ the foliation JF^. has a singular point of saddle type, to which we associate 
a Dulac map (or the transition map), as on fig|2] (i). More precisely, for all 
sufficiently small \e\ the fohation JF^ has two separatrix solutions, which are 
transversal analytic curves, depending analytically on e. We may suppose 
that they are the axes {x = 0} and {y = 0} as on fig|2] (n). Let a, r be two 
complex cross-sections (complex discs) to the two separatrices, parameterized 
hy z = f\a and z = f\r- In these coordinates the Dulac map is the germ of 
analytic map 

4:(M+ 0)^(M+,0) 

defined as follows: ii z e crR IR+ then ds{z) E rCiRf is the intersection with r 
of the orbit 7e(^) of ([2]), passing through z E a. This geometric definition of 
de allows to control to a certain extent its analytic continuation in a complex 
domain. 

3.1 Analytic continuation 

The Dulac map is analytic and hence allows an analytic continuation on 
some open subset of the universal covering a, of a\ {0}. The domain of the 
continuation depends on e, and obviously d^i^z) = z. 

Let us parameterize the universal covering a, by polar coordinates p > 

0,(p EM., z = pexp ip. 

Theorem 2 There exists Eq > and a continuous function 

if ^ p(y?) 
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such that the Dulac map allows an analytic continuation in the domain 

{e, p,yD) e C X a, : \e\ < eo,0 < p < p{<p)} (6) 

The proof of the above Theorem in the 0-parameter case is well known, and 
in the multi-parameter case it is the same. For convenience of the reader 
it will be given in Appendix |^ This proof shows even more: the analytic 
continuation of the Dulac map in the domain ^ can be accomplished in a 
geometric way as follows. 

Let {7o(-2)}j;, 70(2;) C be a continuous family of loops connecting a, 
to r,. For z e a (1 IR+ we suppose that 70 (-2) is the real orbit oi df = 
contained in the first quadrant x > 0,y > 0. We note that, although the 
family {■jo{z)}z is not unique, the relative homotopy class of each loop 70 (-2) 
is uniquely defined for all z E a,. It follows from the proof of Theorem [2] 
that {7o(-2)}z allows a deformation to a family of paths {7e(z)}2, connecting 
0", to r,, tangent to the leaves of JF^, and defined for all e, z in the domain 



3.2 The Poincare-Pontryagin integral 

In what follows a crucial role will be played by the integral f^^^ uq, and its 
generahzations. Namely, let 

K C {(p,^) G a. : < p < p{^)} 

be a compact set, where p{ip) is as in Theorem 2 As ^.^^ Ue is continuous 
in z, e and 



ds{z) - z = j df = e uje 

then the following Lemma holds 
Lemma 1 (Pontryagin|16|) 

d,{z) = z + e! ujo + 0{e^) (7) 
-^70(2) 

uniformly in z E K. 
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The function ujq is the so called Poincare-Pontryagin integral associated 
to the deformed foliation JF^. It follows from the argument principle that, 
that if |e| is sufficiently small, the number of the zeros of de{z) — z in the 
compact K is bounded by the number of the zeros of the Poincare-Pontryagin 
integral /^j^^^) in K (counted with multiplicity). It might happen, however, 
that the Poincare-Pontryagin integral vanishes identically. In all cases there 
is an integer d > 1 and an analytic function Md 7^ in a neighborhood of K, 
such that 

de{z)=z + e''M,iz) + 0i6''+') (8) 

uniformly in z E K, provided that the Dulac map is not the identity map. 
is the so called higher order Poincare-Pontryagin function and its zeros in K 
bound as before the number of the zeros of ds{z)—z. As we already mentioned 
in section [2} there is an integral representation for Md as an iterated integral 
of length at most d along 70(2;). 

Our aim is to obtain a bound for the zeros of de{z) — z in a domain 
K which is open and connected. Even if the estimate ([T]), ([s]) do allow an 
extension to such a domain K, the argument principle can not be directly 
used. For this purpose we consider rather the imaginary part of the Dulac 
map. 

3.3 The zero locus of the imaginary part of the Dulac 
map 

We shall describe the zero locus of the imaginary part of the Dulac map d^ 
in an appropriate sector 

ne = {zeC:lmd,{z) = 0}r\V (9) 

P = {0<p<p(<^),0<<^<y}. 

The surprising fact about He is that it is a smooth real-analytic plane curve 
in P C = C. Even better, the curve He can be conveniently approximated 
in terms of higher order Poincare-Pontryagin functions. 

The foliation Fq has a first integral defining a fibration with fibers Ft, see 
([1]). Let 6{t) C Ft be a continuous family of closed loops 6{t) C Ft vanishing 
at the saddle point when t tends to 0. The orientation of 6{t) is chosen as 
follows. Let 7o(t) be the family of loops defined in the Appendix [Aj For real 
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positive t they coincide with real orbits of J-'o connecting o""*" to r"*". Then, 
the homotopy classes of 70, 5 satisfy 

7o{ten-io{te~n = m- m 

Therefore, the exact orientation of 6(t) can be computed by the Picard- 
Lefschetz formula (but we do not need this). 

Let r be, as before, a cross section to the fiber Fq, see fig. [2j Consider 
the holonomy map /i| associated to the family {S{t)}t and to the deformed 
foliation JF^ 

hi -.T T 

.d 



z ^ z + e'^MJz) + . . . 



The anti-holomorphic involution 

{x,y) ^ {x,y) 
induces, for t G M an anti-holomorphic involution 

Ft^Ft 

which on its turn sends the free homotopy class of the loop 6{t) to the class 
of —6{t). Therefore the function Ma is pure imaginary for real values of t. 

Lemma 2 The zero locus He of the imaginary part of the Dulac map is a 
smooth real-analytic curve o/M^ = C of the form 

= {z = u + iv:v = —Md{u) + e), n < 0} n P (12) 

Zi t 

where R{u, e) is an analytic function. 

The above Lemma is the main technical result of the present paper. The 
analyticity of the zero locus He is responsible for the algebraic-like behavior 
of the Dulac map. 

Proof of Lemma [2j Consider the cross-sections (complex discs transversal 
to the separatrices) cr^, as shown on figjsj simultaneously parameterized 
as before by the restriction z of the first integral f{x,y) on them. The cross- 
sections a,T shown on fig|2]are denoted, from now on, by a~^,T~^. Denote 

^>o = n {(x, y) : fix, y) > 0}, a+Q = a+ n {{x, y) : /(x, y) < 0} etc. 
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Figure 3: The zero locus of the imaginary part of the Dulac map 

Let {7e(-2)}2 be the continuous family of paths, defined in Appendix |A] The 
point 2 G a"*" belongs to the zero locus He if and only if the end of the path 
7e(2;) belongs to r^g- Therefore, such a path allows a decomposition in a 
product 

7e(z) = ae{z) Opeiz) 

where ae{z) is a path connecting 2; G a"*" to a point on ct^q and Pe{z) is a 
path connecting the latter point to a point on t^q, where P^iz) C M^. It 
follows that He is the image of cr^g under the holonomy map 

This already proves that the closure of He C is a smooth analytic curve. 
Once having said this, it is clear that He can be conveniently parameterized, 
which we do next. 



To prove (12), let us note that for t G M the two ends of the loop Q;o(t) are 



real and hence this also holds true for the complex-conjugate loop ao(t). The 
loop ao{t) o aQ^{t) is therefore closed and is homotopic to 5{t) C Ft defined 
above. 
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More explicitely 



2iv 



(hl^, o hl^ - id) o it) 



(/l^-irf)o/ii^(t). 



This, together with (11) and 



h^{t) = t + 0(e), = t + 0(e) 



imphes 



u 



t + 0{e) 



2i 



(Mrf(t) + 0(£)) 



where, by abuse of notation, 0{e) means a function analytic in t,e, which 
vanishes identically for e = 0. This proves the identity (|9]).A 



4 Cyclicity of two-saddle cycles 

In this section we prove TheoremjTJ Using the notations introduced in section 
[2} we suppose that the vector field Xq has a two-saddle loop T2 and an 
analytic first integral / in a neighborhood of it. We suppose that / has 
Morse critical points at the two saddle points pi, p2 of Xq. Consider the 
Dulac maps dl, d^ associated to the corresponding fohation, see figjij Each 
map d\ is a composition of a "local" Dulac map (as in section [s]) and two 
holomorphic holonomy maps. From this it follows that Lemma [2] applies to 
dl, i = 1, 2, too. We parameterize each cross-section by the restriction z of 
/ on it. The function d\, i = 1,2, is multivalued and has a critical point 
at Si{e) G M, Si{0) = 0. The functions Si are real analytic. We consider 
first the case e > and we may suppose that Sj(0) = 0, Si{e) < S2{e) for all 
sufficiently small e, (the case e < is studied in the same way). Our aim is to 
bound the number of those zeros of the displacement map dl — d"^ which are 
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Figure 4: The Dulac maps dl and 



real, bigger than S2{e) and tend to as e tends to 0. Note that these zeros 
correspond to the fixed points of the Poincare first return map P^. Indeed, 

dl — dl = dl o ((Pg — id), where dl = id. (13) 

We shall count the zeros of the displacement map in the larger complex 
domain of the universal covering of C \ defined as follows. It is 

bounded by the circle 

Sr = {z: \z\ = R}, (14) 

by the interval [si(e), S2(e)], and by the zero locus of the imaginary part of 
the Dulac map dl for 3?(z) < si{e), as it is shown on figjs} The numbers 
e, R are subject to certain conditions explained bellow. The zeros of an 
analytic function in a complex domain equal the increase of the argument of 
the function along the border of the domain, divided by 27i (the argument 
principle). To bound the increase of the argument we shall count the number 
of the zeros of the imaginary part of the function, along the border of the 
domain. 

Choose first the real numbers eo,R> as follows. Let 

dliz) - dliz) = e'^Mdiz) + 0(e'^+i^ 



and let z^ilogzY be the leading term of Md{z). Then, by (13), v = z^(Pe) 
is the characteristic number of the Poincare map P^ associated to T2. We 
choose P > so small, that the increase of the argument of Md{z) along the 
circle Sr is sufficiently close to the increase of the argument of z^{\ogz)n 
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Figure 5: The domain 



along Sr. We fix R and choose eo > so small with respect to R, that for all 
e, |£:| < Eq^ the increase of the argument of dl{z) — dl{z) along the circle Sr 
is sufficiently close to the increase of the argument of Md{z) along Sr. This 
is indeed possible, according to Lemma [T) The conditions that we impose on 
Sq^s and R will be denoted (by abuse of notations) as follows 

1 >> R » eo> e > 0. 

To evaluate the increase of the argument of dl{z) — dl{z) along the interval 
[si{e), S2{e)], we bound the zeros of its imaginary part which equals (along 
the interval [si(e), S2{e)]) to the imaginary part of —d'^{z). In other words, we 
need to estimate the number of intersection points (counted with multiplic- 
ity) between the zero locus of the imaginary part of d'^{z) and [si{e), S2{e)]. 
According to Lemma [2] this number of intersection points is bounded by the 
multiplicity of the zero at the origin of the Poincare-Pontryagin function of 
the holonomy map h^^. This multiplicity equals i^{hgj. 

Finally, we arrive at the most delicate point in the proof of Theorem [T] : 
evaluate the increase of the argument of dl{z) — d'^{z) along the zero locus of 
the imaginary part of the Dulac map d] for 3?(z) < si{e). For this purpose 
we bound the zeros of the imaginary part Im {dl{z) — d'^{z)), along the zero 
locus of lmdl{z). Thus, we need to estimate the number of intersection 
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points (counted with multiplicity) between the zero locus of the imaginary 
part of d1{z)^ and the zero locus of the imaginary part of d\{z). 

Recall that to a Dulac map we associated a family of vanishing loops 



5i{z) with orientation prescribed by (10). With this convention, the orien 



tation of the loops 5i and 82 do not agree : if ^{t) is the family of periodic 
orbits of jFg for t > 0, then the intersection indices of the homotopy classes 
of 5i{t),52{t) with 7(t) have opposite signs. In order to have the same con- 
vention as in the formulation of Theorem [T] we reverse the orientation of 82- 
With this convention if 

hl^:z ^ z + e'^^ Mi{z) + . . . (15) 

then, by Lemma [2| the zero locus of the imaginary part of the holonomy 
maps hl^ , is given by 

{z = u + w:v = —Mdi{u)+e'^'+^R\u,€),u< 0}nV (16) 

{z = u + iv:v = -—Md2{u) + e'^^+^R^{u,e),u < 0}nV (17) 
^ c 

respectively, where R^,R^ are appropriate analytic functions. 

We conclude that the number of intersection points of the above analytic 
curves coincides with the multiplicity of the zero at the origin of either M^i 
(if < d^), or Md2 (if d^ < d^), 01 M^i + M^z (if d^ = d^). The number of 
intersection points equals therefore to the characteristic number z/(/i|^ ° ^la)' 

Summing up the above information we get that the increase of the argu- 
ment of dl{z) — d1{z) along the boundary of the complex domain is not 
bigger than 

v{P,) + v{hl) + v{hlohl^) + 2. 

The above estimate can be slightly improved, by taking into consideration 
the fact that the imaginary part of d\{z) — dl{z) vanishes at S2{e). Theorem 
[T]is proved. □ 



A Proof of Theorem 2 



In this appendix we will prove Theorem [2]in the slightly more general context 
of multi-parameter analytic deformations. This will be used in Appendix [B| 
Note that the zero-parameter case is well known [T3l 119]. 
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Figure 6: The path l{z) 



Consider a A^-parameter analytic family of analytic plane vector fields 
Xx, A G (]R^,0), such that Xq has a hyperbolic singular point at the origin. 
It is known, since Briot and Bouquet, that Xq has two transversal invariant 
analytic curves which can be supposed to coincide with the = and 

y = 0, that is to say 

Xo = Aix(l + ...)— + \2yil + ...)^, A1A2 < 0. 

The proof is as follows : a formal change of the variables first removes some 
(but not all) non-resonant terms, and then one verifies the convergency of 
the transformation, see [HI Appendice II]. Exactly the same proof apphes, 
however, to the family X^. One can show in this way that Xx is analytically 
orbitally equivalent to the following (slightly improved) normal form 

d , . d 

where r = r(A), a = a{x,y,X) are appropriate analytic functions in their 
arguments, r(0) < 0, see [HI Appendice 1]. We shall suppose, without loss 
of generality, that there exists a constant c > such that r, a are analytic in 
the complex domain 

= {x, A) : < 2, |?/| < 2, |A| < c}. 
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After a further linear re-scaling of x,y, we may suppose that \a{x,y, A)| is so 
small in Vc, that 

r(A) + xy.a{x, y, X) ^ . 
After this preparation, choose the cross-sections 

a = {y = 1},T = {x = 1} 

and consider the corresponding Dulac map 

d\ : a T 

Z dyXz) . 

To prove Theorem [2] we have to show that the constant c > can be chosen 
in such a way, that for every v^o > there exists < < 1, such that the 
Dulac map allows an analytic continuation in the sector 

{z G C : |2;| < zq, \arg(z)\ < Lpo} x {A : |A| < c}. 

The proof is similar to the proof of [T9l Theorem 7], the only difference being 
the presence of the parameter A. We shall construct a continuous family of 
paths 7a (2;) contained in the leaves of the foliation J^x defined by the vector 
field Xx in C^. Each path 7a(-2) starts at the point {x = z^y = 1) and ends 
at the point {x = l,y = dx{z). The path 7a(-2) is constructed by lifting the 
path l{z) contained in the x-plane {y = 0} and shown on figjoj with respect 
to the projection 

7r:C2 ^ C 
{x,y) H-^ X 

Indeed, the fohation J^x is transversal to the projection vr except along 
the leaf x = 0, provided that r + xy.a 7^ 0. The resulting path path ■jx{z) is 
shown on fig. [7| 

To prove the existence of 7a(2;), consider the solution y = y{x) associated 
to the vector field Xx, with initial condition y{z) = 1. We have to show that 
the solution y = y{x) exists when x is restricted to the path l{z). The path 
l{z) is composed by an arc and a segment. We consider them separately 

• Along the arc 

X = l^le**^, y = 0,0 < (p < arg{z) 
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V 




parameterized by if we have 

dy = —y{r + \z\e^'^y.a)dip 
d\y\ = —\y\.\z\lm.{e'''''^y.a)dip. 

Therefore, if |z| < 2:0 is sufficiently small, then |?/(|z|e*'^)| < 2 when 
< < arg{z) 

• Along the segment 

X E [\z\,l] 

we have similarly 

xdy = —y{r + xy.a)dx 
xd\y\ = —\y\{r + ^{xy.a))dx. 

The derivative ^ is therefore negative, the function \y\{x) decreasing, 
so \y\{x) < 1. 

□ 

B Mult i- parameter deformations of Hamilto- 
nian two-saddle loops 

Consider, as in the preceding Appendix, a A^-parameter analytic family of 
analytic plane vector fields Xx, A G (]R^,0). We suppose that Xq has a 
Hamiltonian two-saddle loop r2 bounding a period annulus. This case is 
easily reduced to the one-parameter case studied in the present paper by 
making use of a standard procedure based on the Hironaka desingularization 
theorem. In this Appendix we indicate the main steps. 

B.l Principalization of the Bautin ideal 

In this section we follow [8|l20]. Let zq E a, zq ^ T2 and consider the Poincare 
map 

00 

Px{z) = z + Y,ci^wi^-^oy■ 
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The Bautin ideal, associated to Pa, is the ideal J =< aj >, generated by 
the germs of the analytic coefficients ai{.). It is Noetherian, so generated by 
a finite number of coefficients, and moreover does not depend on the choice 
of Zo [19]. More generally, let {7(^)} be any continuous family of closed 
loops in the fibers of the fohation {Fq, intersecting the cross-section a. For 
all sufficiently small A the (germ of) holonomy map 



: a a 



is defined. In the same way, we associate to h^{z) a Bautin ideal, generated 
by the coefficients aj(A) of the expansion of h^{z) with respect to 2; — zq. As 
before it is Noetherian, and does not depend on the choice of zq (with the 
same proof). 

We may assume, without any loss of generality, that the Bautin ideal is 
principal. For this we use a variant of Hironaka's desingularization theorem 
as follows. 

Let ipo, ipi, . . . ,ipphe non-zero analytic functions on a smooth complex or 
real analytic variety X. The indeterminacy points of the rational map 

(f-.X—^ pp 
can be eliminated as follows [HI [2] 

Theorem 3 (Hironaka desingularization) There exists a smooth ana- 
lytic variety X and a proper analytic map n : X ^ X such that the induced 
map (f = (f o TT is analytic. 

X 

X - - ^ pp 



Let Ox be the sheaf of analytic functions on X and consider the ideal sheaf 
/ C Ox generated by ipo, ipi, . . . , ipp. The inverse image ideal sheaf of I under 
the map vr : X — > X will be denoted vr*/. This is the ideal sheaf generated 
by the pull-backs of local sections of I. We note that vr*/ may differ from 
the usual sheaf-theoretic pull-back, also commonly denoted by vr*/. A simple 
consequence of Theorem [3] is the following 

Corollary 1 The inverse image ideal sheaf n*I is principal. 
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This is called the principalization of I. Indeed, as the induced map ip is 
analytic, then for every A G X there exists j, such that the functions (fi/^pj, 
i = 1,2, are analytic in a neighborhood of A. Therefore there is a 

neighborhood U of X such that divides in the ring of sections Ojy 
of the sheaf (9^^, that is to say is generated by (pjljj. 

In our context X = a is the cross-section to the family of periodic orbits 
{■y{z)} and ipo, v^i, ■ ■ ■ ,fp are the germs of analytic functions which generate 
the Bautin ideal associated to the holonomy map h^. To apply Theorem [s] 
we assume that a is a polydisc on which ipi are analytic, and the divisors 
(ipi) intersect transversally the boundary of a. After applying the Hironaka's 
theorem, the origin of cr is replaced by a compact divisor 7r^^(0), along 
which the inverse image ideal sheaf vr*/ is principal. 

Suppose now that we have a holonomy map associated to another 
family of periodic orbits 6{t). Let ipo, ipi, . . . ,'ipq be generators of the corre- 
sponding Bautin ideal J. As before we assume that ijji are analytic on a, with 
divisors transversal to its boundary. Applying twice the Hironaka's theorem 
we get a new analytic variety smooth analytic variety X and a proper ana- 
lytic map n : X ^ X such that the induced maps = (y9 o vr and ip = ifj o n 



are analytic, see diagram (19). 




(19) 



X - - ^ pp 



The inverse image ideal sheaf vr*/ and ti* J are both principal along the 
compact divisor 7r~^(0). 

B.2 Multi-parameter version of Theorem [l] 

We begin by formulating the multi-parameter version of Lemma |2| let {5{t)] 



be the family of vanishing loops defined in section 3^ In agreement with 
the preceding section, let us suppose that the Bautin ideal of the holonomy 
map is principal. We have, therefore (compare to (11) 



^'^■■^ ^ ^ . . (20) 

z ^ z + (^(A)(M(z) + i?(z,A)) ^ > 

where R{., .) is analytic, R{z^ 0) = 0, and M(.) is the highest order Poincare- 
Pontryagin function. 
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Px = h'hl,hlh'ohi (22) 



Lemma 3 The zero locus He of the imaginary part of the Dulac map is a 
smooth real-analytic curve o/M^ = C of the form 

T-l^ = = u + iv:v = ^^{M{u) + R{u, X)),u<0}nV (21) 
where R{u, A) is an analytic function, R{u, 0) = 0. 

The proof of the above Lemma is completely analogous to that of Lemma [2] 
and is therefore omitted. 

In the proof of Theorem [T] we used four Bautin ideals associated to the 
holonomy maps 

D — ,, ,, ,, , , ,, 

''7' "'Si ' ""(52' Si "'S2 

where {7(2;)} is the family of periodic orbits associated to the annulus, 
{5i(z)}, {52(2)} are the vanishing famihes of loops associated to the sad- 
dle points. After an appropriate blow up tt we may suppose that the inverse 
image ideal sheafs of the corresponding four Bau tin ideals are principal along 
the compact divisor 7r~^(0) C X, see section [bT Let A be a local variable on 
the smooth variety X. The cychcity Cycl{T2, {^x, ^x^)) is the maximal num- 
ber of limit cycles which JF^ can have in an arbitrarily small neighborhood 
of when A tends to Aq- Denote also 

Cyd{r2,J'x) = Cycl{r2,{J'x,J'o). 

Clearly 

Cyd{T2,Tx)= sup Cyd{T2,iJ^x,J^xJ) 

A67r-i(0) 

and because of the compactness of 7r"^(0), there exists Aq G vr"^(0) such that 
Cyd{T2,J'x) = Cyd{T2, {:F-„J^-J). 

The above considerations show that, without any harm, we may sup- 
pose that X = cr, A = A, Aq = 0, and the Bautin ideals associated to the 
holonomies (22) are principal. Consider the circle Sr defined in (14). The 
Bautin ideal of the Poincare map Pa coincides with the Bautin ideal of the 
displacement map d\ — d\ and 

dl{z)-dl{z) = ^{X){M{z) + R{z,X)) 

where cp is the generator of the Bautin ideal, R is analytic, R{z, 0) = 0, and 
M{z) is the Poincare-Pontryagin function. As before M satisfies a Fuchs 
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equation with a singularity at z = 0. We choose R so small, that the increase 
of the argument of M{z) along Sr, arg(z) < tt is close to the increase of the 
argument of the leading term of M. We note that if z'^ {\og{z))^ is the leading 
term of M, then 

where e ^ \{e) is a one-parameter deformation (a specialization), such that 
(p{X{.)) 7^ 0. We fix R and choose > so small, that for all A, such that 
|A| < Eo, the increase of the argument of the displacement map d\{z) — d\{z) 
along Sr^ arg(z) < vr, is close to the increase of the argument of M{z). By 
making use of Lemma |3} the proof of Theorem [T] is completed as in section |4] 
Note that the characteristic numbers which appear in the estimate of the 
cyclicity may be obtained as characteristic numbers corresponding to a one- 
parameter analytic deformation J-'\(e) of J-'o 

e \{e), A(0) = 0, lim A(£) = Aq 

provided that the generators of the Bautin ideals do not vanish identically 
along this deformation. Therefore, a multi-parameter version of Theorem [T] 
can be formulated as follows 

Theorem 4 There exists a germ of analytic curve e — A(e), A(0) = m the 
parameter space, such that 

CycliV^, ^a) < 1 + ^{Pxie)) + max{i.(/if )), v{h]^^)} + v{hf^ o /.Jf ). (23) 
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